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Abstract 

Parity-odd asymmetries in the electromagnetic decays of compound nuclei can sometimes be 
amplified above values expected from simple dimensional estimates by the complexity of compound 
nuclear states. In this work we use a statistical approach to estimate the root mean square (RMS) 
of the distribution of expected parity-odd correlations • k™, where is the neutron spin and 
& 7 is the momentum of the gamma, in the integrated gamma spectrum from the capture of cold 
polarized neutrons on Al, Cu, and In and we present measurements of the asymmetries in these 
and other nuclei. Based on our calculations, large enhancements of asymmetries were not predicted 
for the studied nuclei and the statistical estimates are consistent with our measured upper bounds 
on the asymmetries. 
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I. INTRODUCTION 



One might assume that a quantitative treatment of symmetry breaking in neutron reac- 
tions with heavy nuclei would not be feasible. However, theoretical approaches exist which 
exploit the large number of essentially unknown coefficients in the Fock space expansion of 
complicated compound nuclear states in heavy nuclei to perform calculations that can be 
compared to experiment. If we assume that it is possible to treat the Fock space components 
of the states as independent random variables, one can devise statistical techniques to calcu- 
late, not the value of a particular observable, but the root mean square of the distribution of 
expected values. This strategy has been used successfully to understand certain global fea- 
tures of nuclear structure and reactions Q] . The distribution of energy spacings and neutron 
resonance widths, for example, has been known for a long time to obey a Porter-Thomas 
distribution |2| in agreement with the predictions of random matrix theory, and statistical 
approaches have been used to understand isospin violation in heavy nuclei [3|. 

The complexity of the compound nuclear states can also amplify the size of the parity- 
odd asymmetries by several orders of magnitude relative to single-particle estimates. This 
large amplification makes it practical to use nuclear parity violation, generically expected on 
dimensional grounds to possess amplitudes seven orders of magnitude smaller than strong 
interaction amplitudes, as a new setting to investigate the validity of these statistically-based 
theoretical approaches. Statistical analyses have successfully been applied recently to an 
extensive series of measurements of the parity-odd correlation s* n -p n in the A=100-200 mass 
region in neutron-nucleus scattering performed at Dubna, KEK, and LANSCE 
Although the comparison between theory and experiment in this work is still hampered 
somewhat by the lack of precise knowledge of the weak NN amplitudes and their possible 
modifications in the nuclear medium, theory and experiment appear to be in agreement 
at about the 50% level. Given the extreme complexity of the states involved, agreement 
between theory and experiment at this level must be counted as an overall success for the 
statistical approach. 

Parity violation in the gamma decays of nuclei is another example where statistical meth- 
ods may be employed to estimate observables. In this case the observables involve the 
parity-odd correlation s* n ■ fc 7 where s* n is the neutron spin and fc 7 is the momentum of the 
gamma 0,0,0]. Just as for neutron scattering, neutron capture on elements with a large 



3 



number of nucleons produces compound nuclei in highly excited states. These nuclei exhibit 
a huge (> 10 6 ) number of possible state configurations with different angular momenta and 
parity and the number of transitions with different amplitudes that the compound nucleus 
may make to its ground state is correspondingly large as well. Because of the large number 
of energy levels in the compound nucleus, formed by neutron capture, one may hope that 
the calculation of the mean square matrix elements for the transition amplitudes may also 
amount to a summation of a large number of uncorrelated random contributions as in the 
case of the total cross section. One can then use statistical arguments to estimate the RMS 
value of the parity-odd 7-ray asymmetry. 

However the case of parity violation in (11,7) reactions in heavy nuclei is not quite as 
simple as parity violation in the total cross section for both theoretical and experimental 
reasons. For the total cross section, the amplification of parity violation effects is dominated 
by the mixing amplitude of the weak interaction, between two compound nuclear states 
of opposite parity (in practice s-wave and p-wave compound states). Since the total cross 
section is proportional to the forward elastic scattering amplitude, by the optical theorem, 
there is only one such contribution for any pair of opposite-parity compound states. For 
inelastic processes such as the (11,7) reaction, however, the weak mixing between compound 
states can occur in either the initial or final nuclear states, and since these states are distinct 
in an inelastic reaction there are two possible sources of compound nuclear amplification of 
the parity-odd effect rather than one ^. Because of the large density of states in the initial 
state near neutron separation energy, the initial state mixing will involve a larger number 
of components in the wave function for gamma transitions to low-lying states and therefore 
lead to a larger amplification. However one also has contributions from transitions to higher- 
lying states where final-state mixing is somewhat more important. Experimentally, precise 
measurements of parity-odd asymmetries are more practical for the total integrated gamma 
spectrum rather than individual gamma transitions. But a calculation of the asymmetry 
of the integral gamma spectrum requires an additional averaging over the large number 
of distinct final states. In addition the integral measurement also senses gamma cascades 
in addition to single transitions. Parity-odd correlations in the integrated gamma spectra 
of 35 C1, 81 Br, 113 Cd, 117 Sn, and 139 La have previously been calculated by Flambaum and 
Sushkov (| and by Bunakov et. al. j|. However, more experimental information on parity- 
odd asymmetries in integral gamma spectra from heavy nuclei are needed in any attempt 



4 



to make progress in this area. 

We have searched for parity-odd directional 7-ray asymmetries in the capture of cold 
polarized neutrons on 27 Al , Cu, and 115 In at the Los Alamos Neutron Science Center 
(LANSCE). We have performed a simple statistical estimate of the mean square value for 
the parity odd asymmetries in these nuclei and obtain expected upper bounds which are 
consistent with experiment. In addition, we performed measurements of the directional 7- 
ray asymmetry for polarized cold neutron capture on 35 C1 and on 10 B. 35 C1 is known to 
possess a large parity-odd gamma asymmetry and it is used to verify the sensitivity 

of our apparatus. 10 B is used extensively throughout the experiment, for neutron shielding. 
Searches for parity-odd gamma asymmetries on several other nuclei are in progress. 

These measurements are being conducted in preparation for an experiment to search for 
the parity violating 7-ray asymmetry in the capture of polarized neutrons on protons by the 
NPDGamma collaboration. The apparatus constructed for this measurement is capable of 
measuring 7-ray asymmetries with an accuracy of 10~ 8 . 

The remainder of the paper is organized as follows. We first provide a brief theory section 
in which we outline the calculation and estimate the expected root mean square of the 7- 
ray asymmetry in a current mode 7-ray detector from the nuclei used in the experiment. 
We give a short overview of the experimental layout and then describe the measurements. 
We conclude with a discussion about the results and the associated implications. 



II. THEORY AND STATISTICAL ESTIMATES 



The simplest nuclear reaction which can produce a parity-odd directional distribution 
of 7-rays is the capture of polarized neutrons on protons. The differential cross-section 
in this simple system can be calculated explicitly from the transition amplitudes of the 
electro-magnetic part of the Hamiltonian between initial (capture) and final (bound) two 
nucleon states, which possess mixed parity due to the NN weak interaction. In the n + p — > 
d + 7 reaction the primary process is the strong interaction induced parity conserving Ml 
transition between the singlet and triplet S-wave states: 1 Sq , 3 Si. The weak interaction 
introduces a small parity non- conserving admixture of P-wave states in the initial singlet 
and the final triplet S-wave states. The largest contribution to the hadronic weak interaction 
comes from pion exchange and the measurement of the parity- violating up-down asymmetry, 
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Ay, in the angular distribution of 2.2 MeV 7-rays with respect to the neutron spin direction 
(Eq. [IJ , almost completely isolates the term proportional to the weak pion-nucleon coupling 
constant /„. jl^ . 



— oc — (1 + A y cos9) (1) 



da 1 

Here cos 9 is angle between the neutron spin direction and the 7-ray momentum. 

For the n + p — > d + 7 reaction, it can be shown that there is a simple expression for the 
7-ray asymmetry in terms of the matrix elements between initial and final states 

A o:Rc e{3PllEll3Sl) (2) 



Here 



(3) 



AE 

and a = {J, L, S,p} (p = parity). 

In heavy nuclei the interference term which produces the asymmetry is much more com- 
plicated, involving many states. Here, a neutron may capture into an S or P wave state close 
to the neutron separation energy (S n ) and the weak interaction mixes the corresponding am- 
plitudes perturbatively. For almost all nuclei except in few body systems it is essentially 
impossible to calculate the parity violating asymmetry from the strong and weak Hamil- 
tonian, because of the large number of 7-ray transitions. However, because of the large 
number of possible electromagnetic transitions in the compound nucleus the calculation of 
the mean square matrix elements for the transition amplitudes amounts to a summation of 
a large number of uncorrelated random amplitudes which are approximately independent of 
the transition energy for E < S n . One can then hope to use statistical arguments to estimate 
the RMS value of the asymmetry from nuclei close to a certain neutron separation energy. 

Due to the large density of states close to the neutron separation energy and the corre- 
spondingly small level spacing D ~ AE C , parity violation is expected to be dominated by 
the mixing of the two closest S and P wave states near S n , in the initial or capture state, 
and it is expected on general grounds that parity violation due to mixing with lower lying 
states may be neglected. The parity violating asymmetry comes from interference between 
El and Ml gamma transitions. The 7-ray asymmetry from the decaying compound nucleus 
as measured in a current-mode gamma detector is given by 
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'j2 Jf (J}\El\J?)(Jf\ML\J$E*v 






[\(J* f \Ml\J?)\> + \(Jf\El\jf)\A 





A, = e • £ • F(J T , Jf) / , x J • (4) 



Here the transitions are between initial (i) and final (f) compound nuclear states with total 
angular momentum (Jj, Jf) and parity (p,p 7 ). F(Jt, JiHs the angular momentum coupling 
factor resulting from the compound state polarization 6[ : 



F(J T ,J i ) = (-l) 2J ^ 2+Jr 3(2J i + l) 



1 1/2 1/2 

Jt Ji Ji 



Jt is the angular momentum of the target nucleus before neutron capture. 

The dependence on the 7-ray transition energy E^ in eqn. (0J) comes from the phase-space 
factor (E^ 2 ) in the transition amplitude and the linearity (oc E-y) of the detector response 
as a function of energy in a current mode gamma detector. The factor 

_ J2fI-Y,if E i,if _^ 1 Jo" E^ l p f {E 1 )dE 1 



arises because the current mode gamma detector possesses no energy resolution and therefore 
sees a superposition of currents from all transitions. This has the effect of diluting the 
asymmetry (0 < £ < 1). Here, 



(K^|Ml|Jf)| 2 + |(^|El|Jf)| 2 )^. 



E./K^IMiijnp + KjjiEiiJf)! 2 )^ 

is the relative intensity of a given transition. 

We estimate the density of final states using the Back-Shifted Fermi Gas model (BSFGM) 



as 
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14] 



( F \ _ \T 2 J + 1 exp[2y/a(£ a - A) - J(J + l)/2a 2 } 
PAx) ZfUyftoZaV* (E x - A + t)5/4 ' W 

where J is summed over Jf — 1, J/, Jf + 1 for each final compound nuclear state. Here, 
a [MeV -1 ] and A [MeV] are determined from experimental data and the temperature pa- 
rameter t is defined by E x — A = at 2 — t. a 2 = I e fft/h 2 ~ 0.015A 5//3 t is the spin cut-off 
parameter and the effective moment of inertia I e ff takes on values between 50% and 100% 
of the rigid body moment of inertia I r i g = ^MR 2 . The level density (Eq. is derived 
assuming random coupling of angular momenta and the spin cut-off parameter arises as a 



result of this treatment [14j . The excitation energy E x = S n — E 1 is the energy of the nucleus 
after the 7-ray transition from the capture state. E x may be zero if the transition is to the 
ground state. Figure [T] shows the predicted density of final states, using this model, for the 
27 Al , Cu, and 115 In, nuclei as a function of 7-ray energy. 
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FIG. 1: Density of final states in the excited compound nuclei investigated in this work as a 
function of 7-ray transition energy up to the neutron separation energy. There are many more 
states at low gamma energies than at high energies, and the decaying nucleus emits many low 
energy gamma-rays before reaching the ground state. The level density is calculated according to 
the Back-Shifted Fermi Gas Model. 

For comparison and to estimate the uncertainty in the calculated asymmetries due to the 
model we also determine the asymmetries using a constant temperature model (CTM) of 
the final density of states [3] 

Pf (E x ) = ^exp(E x -A)/T . (6) 

The aim of this calculation is to find a simple "generic" formula that holds for many 
nuclei and provides a good estimate of the size of an asymmetry one can expect in a mea- 
surement of this nature. The denominator in eq. 0] is the parity allowed transition from 
the initial compound state ( Jj) after capture of an S-wave neutron. This transition has the 
largest amplitude and basically determines the intensity of the gamma signal. In general El 
transitions outnumber Ml transitions and the denominator is primarily El for most nuclei. 
We point out though, that if one has initial states such that all or most parity allowed 
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transitions are Ml in the range of expected 7-ray energies, as determined by the density of 
states (as is the case for Al and In), then the denominator would be Ml. 

The root mean square of the 7-ray asymmetry can then be estimated as follows: We use 
the electric and magnetic dipole transition rates which are given by 



El 



2tt 



(Jf I El I Jf) ) Pf (S n ) 



T M1 = 27r(\(J p f \Ml\Jf}n Pf (S n ) 



(7) 



respectively. The transition rates are strength functions. As the density of states increases, 
the average matrix element squared decreases and the transition rates are constant or slowly 
varying functions of energy. 

The root mean square of the detected intensity of the gammas that depopulate the initial 
state is given by taking the average of the squared denominator in eqn. 0] Then, invoking 
the randomness in the transition amplitudes (under the assumption that the correlation is 
zero, so that the cross terms vanish), we find 

2 



(|(Jj'|El|Jf)| 2 + |<Jj?|Ml|Jf}| 2 ) El 




E(V/'i E1 i J "> 



y 2np f (S n 



r 2 4- r 2 

1 El + 1 Ml 



S >i 




]T(|<jj|Mi|jf>| 2 W 



S ji 



(8) 



47T 2 p 2 (S n , VJ0 

The factor pj{E 1 )dE 1 arises in the standard fashion, when converting the sum over final 
states into an integral. 

The root mean square of the interference term in the numerator gives 

2 



4^E( J /l E1 l J F)( J fl M1 l J /)^ 

4^(|{j/|El|Jf> 2 \(|(Jf|Ml|Jf)| 2 )< l/ 




El 



El 



Ml 



2-Kp f (S n ) 2irpf(S n 



-Pf(E 1 )dE 1 



~~ 2~72 



7T 



p){S n 



Sn 

E*p f (EjdE^ . (9) 



Where we again used the randomness in the transition amplitudes in going to the second 
line. With this, the mean square asymmetry can be estimated for each target from. 



(M) ~ 2F( J T , J,K 



f Sn E*p f (E 7 )dE, 



\ T Bi + r Mi E* Pf (EjdE,y 



To calculate the root mean square asymmetry for a particular nucleus, one must then 
determine whether the transitions to the ground state are mostly El or Ml and omit the 
corresponding amplitude in the denominator. In the case of 27 Al and 115 In we then have 
S^, while for Cu we have 

J- El 1 Ml 

Substituting the experimental value of the hadronic weak mean square matrix element 
(rv = l-8lg;J x 10 " 7 eV ) B and usin g 

2 _ rv i ^ rv 

6 ~ 27rpi D 2 ~ 2ttD 

together with the fact that El transitions are approximately 10 times faster than Ml transi- 
tions, T^i ~ lOr^i ^ ne ro °^ mean square asymmetry can be calculated for different 
nuclei and neutron separation energies. When evaluating eqn.Elfor aluminum, for example, 
the single particle level spacing is approximately D ~ 120 keV, the dilution factor £ 2 ~ 0.6 
and the ratio of integrals in eqn. ^Jcan be numerically evaluated to give ~ 6.5 x 10~ 2 . The 
expected RMS value of the gamma asymmetry is then about 1.3 x 10~ 7 . The RMS 7-ray 
asymmetry values and other associated variables for the nuclei studied in this work are listed 
111 table □ 



A. Theory Discussion 

The results in tables HI and II VI show no large enhancements. There are several reasons why 
one may expect this behavior. For example, the levels are highly degenerate, the sign of the 
asymmetry is random and the transitions mix incoherently, producing a 1/ suppression. 
There is also no kR enhancement for the direct capture calculations done here, which are 
appropriate for the low neutron energies used in these experiments. 
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In jy] Flambaum and Sushkov calculated the average value for the integral 7-ray spectrum 
relative to the S-wave amplitude for thermal neutrons 
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T2-p(7) 
1 s 1 eff 



4P (E - E s f + \Yl 

which is far from p-wave resonance so that it's contribution to the cross section can be 
neglected. The root mean square asymmetry is given, in their notation, by 



(A 9 ) = -2Re 



F(J T , Ji. 



-T . 



(11) 



E-E p - \iV p ) 3y/2Ji + 1 
Where (r) is an integral over the El and Ml radiative strength functions, detection efficiency 
and density of final states, corresponding to our integral in eqn. |§1 Equation may be 
compared to our result above. Flambaum and Sushkov also state that the 7-ray asymmetry 
arises as a result of the El, Ml interference, that the transitions are random, and that the 
asymmetry is statistically suppressed after averaging. 

The main difference between our calculations and those done by Flambaum and Sushkov 



Calculated RMS 7-ray Asymmetries Values (BSFGM) 
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Calculated RMS 7-ray Asymmetries Values (CTM) 
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TABLE I: RMS 7-ray asymmetry values and associated variables, as estimated from the statis- 
tical pproach (i = f£ n dE 7 E*p(E^)/ ( J^ n dE^E^E^ , D = D Q £ 2Ji + lj , D a was taken 



from 
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is that they consider a p-wave resonance near the thermal (or cold) region mixed by parity 
violation with one S-wave resonance, while our treatment takes account of all S and P-wave 
resonances, but in the tail, far from resonance, at the average spacing D or more. 



III. EXPERIMENT 

The NPDGamma apparatus used for the measurements is located on flight path 12 at the 
Manuel Lujan Jr. Neutron Scattering Center at LANSCE. The LANSCE linear accelerator 
delivers 800 MeV protons to a storage ring, which compresses the beam to 250 ns wide pulses 
at the base. The protons from the storage ring are incident on a split tungsten target at 
a rate of 20 Hz and the resulting spallation neutrons are cooled by and backscattered from 
a cold H 2 moderator with a surface area of 12 x 12 cm 2 . For the measurements described 
here, the cold neutrons were transported to the experimental apparatus by a neutron guide 
and then transversely polarized by transmission through a polarized 3 He cell. Three 3 He 
ion chambers were used to monitor beam intensity and polarization. A radio frequency 
spin flipper was used to reverse the neutron spin direction on a pulse by pulse basis. The 
polarized neutrons then captured on a target placed in the center of the gamma detector 
array. The gamma rays from the neutron capture were detected by an array of 48 CsI(Tl) 
detectors operated in current mode The entire apparatus was in a homogeneous 10 

Gauss field, which was required to maintain the neutron spin downstream of the polarizer, 
with a gradient of less than 1 mG / cm to make spin-dependent Stern-Gerlach steering of the 
polarized neutron beam negligible. 

Figure El shows the flight path and experimental setup. The distance between the mod- 
erator and target is about 22 meters. The flight path 12 beam line consists of a neutron 
guide, a shutter, and a beam chopper. The pulsed spallation neutron source allowed us to 
know the neutron time of flight or energy accurately. The chopper is used to define the time 
of flight frame and to prevent neutrons from different frames to mix and thus dilute the 
neutron energy information. In this experiment the chopper was used to close the beam line 
before the end of the frame which allowed us to take beam-off (pedestal) data for ~ 6 ms at 
the end of each neutron pulse which is needed for detector pedestal and background studies 
(Fig. IHJ) • The last 10 ms after sampling stops is used by the DAQ for data transfer. A de- 
tailed description of the FP12 neutron guide and performance is given in [2^. The measured 
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FIG. 2: Schematic of the experimental setup. 

moderator brightness has a maximum of 1.25 x 10 8 n/ (s ■ cm 2 ■ sr • meV ■ /iA) for neutrons 
with an energy of 3.3 meV. 
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FIG. 3: Normalized signal from the first beam monitor downstream of the guide exit. The solid 
triangles show the signal obtained from a run where the chopper was parked open. The open circles 
corresponds to a run taken with the chopper running. 

The neutrons were polarized by passing through a 12 cm diameter glass cell containing 
polarized 3 He 



ithin). The beam polarization was measured with 
the beam monitors using neutron transmission ( 3 He polarization can be monitored using 
NMR). For gamma asymmetry measurements, the figure of merit is the statistical accuracy 
that can be reached for a certain running time, which is proportional to the product P„,\ /TL. , 
where T n is the neutron transmission through the cell and P n is the neutron polarization 
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The neutron transmission increases with energy whereas the neutron polarization decreases 
with energy. In the analysis of the data the neutron polarization was calculated separately 
for each run by fitting the transmission spectrum to the expression P n = tanh (a c nlPH e ), 
using a 3 He thickness of 4.84 bar • cm, which was separately measured. Here, a c = a a / \[E 
with neutron energy E in units of meV, and a D = 27168 b, nl = 4.84-2.688 x 10 23 atoms/m 2 . 

The primary technique for reducing false asymmetries generated by gain non-uniformities, 
slow efficiency changes and beam fluctuations is frequent neutron spin reversal. This allows 
asymmetry measurements to be made in each spin state for opposing pairs of detectors 
and for consecutive pulses with different spin states, thereby suppressing the sensitivity of 
the measured asymmetry to detector gain differences, drifts, and intensity fluctuations. By 
carefully choosing the sequence of spin reversal, the linear and quadratic components of 
time-dependent detector gain drifts in a sequence can be greatly suppressed. To achieve the 
neutron spin reversal, the experiment employed a radio frequency adiabatic neutron spin 
rotator (RFSR) [2^ which operates at 29-kHz for the 10 G guide field. The neutron spin 
direction is reversed when the RFSR is on and is unaffected when it is off. The spin flip 
efficiency averaged over the beam cross-section ( 5 cm diameter) was measured to be about 
99% 

The polarized neutrons then captured on a target placed in the center of the gamma 
detector array. The targets were thick enough to stop most of the neutron beam by capture 
or scattering with diameters larger than the beam cross section. The capture 7 signals from 
all of the targets measured were large compared to noise and background. 

The housing for the 33 cm 3 liquid CCI4 target was made of Teflon. The CCI4 liquid 
is 99.9% chemically pure, with less than 0.01% water content. The aluminum and copper 
targets consisted of a number of sheets supported by an aluminum frame. Each target sheet 
is an approximately 1 mm thick square with 8.5 cm sides. The arrangement of the target 
into sheets with a gap between the sheets reduced 7 attenuation in the target. The total 
length of the target (including gaps) was 30 cm. Target out background runs and runs with 
the empty frame were conducted as well and the background is taken into account in the 
final determination of the asymmetry (see table |HJ). The boron target consists of a 1 cm 
thick 15 cm by 15 cm sheet of sintered B 4 C glued to an aluminum holder consisting of a 
simple (thin) aluminum sheet. The indium target was approximately 12 mm thick, covering 
a circular cross-sectional area with a radius of ~ 3 cm at the center of the beam. For each 
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Relative Background 



B 4 C 


< 


17% 


Al 


< 


15% 


In 


< 


11% 


CC1 4 


< 


8% 


Cu 


< 


7% 



TABLE II: Targets with their relative background contributions (target in versus target out). In 
each case the maximum value is stated for the detector with the largest background signal. The 
relative amount of background varies, because the magnitude of the 7-ray signal varies with target 
while the target-out background remains constant. 

of the targets the beam was collimated to a diameter of about 5 cm. 





2c in C 
3&tot 


(A* 




Al 


3 x 10' 


-3 


1 


Cu 


2 x 10" 


-2 


0.95 


CCI4 


7 x 10" 


-2 


0.95 


In 


2 x 10" 


-3 


1 


B 4 C 


5 x 10" 


-4 


1 



TABLE III: Spin-flip probability estimate and corresponding corrections to the asymmetry due to 
depolarization. 

The depolarization of neutrons via spin flip scattering from the nuclei dilutes the asym- 
metry. For all targets the neutron depolarization is a small effect which can be estimated to 
sufficient accuracy for nonmagnetic materials using the known neutron coherent and inco- 
herent cross sections. Table ILTT1 lists the estimated spin- flip probabilities for the targets used 
and the corresponding calculated average correction factors (F(ti)). The degree of spin flip 
scattering is neutron energy dependent and a Monte Carlo calculation for the depolarization 
as a function of neutron energy was applied to the data. 

The detector array consists of 48 CsI(Tl) cubes arranged in a cylindrical pattern in 4 rings 
of 12 detectors each around the target area (Fig.HJ). In addition to the conditions set on the 
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detector array by the need to preserve statistical accuracy and suppress systematic effects, 
the array was also designed to satisfy criteria of sufficient spatial and angular resolution, 
high efficiency, and large solid angle coverage [19]. Because of the possible small size of the 
asymmetries and the proposed measurement accuracy the average rate of neutron capture 
and the corresponding gamma rate in the detectors must be high to keep the run-time 
reasonable. Because of the high rates and for a number of other reasons discussed in [lj|, 
the detector array uses current mode gamma detection. Current mode detection is performed 
by converting the scintillation light from CsI(Tl) detectors to current signals using vacuum 
photo diodes (VPD), and the photocurrents are converted to voltages and amplified by 
low-noise solid-state electronics j^. 

In current mode detection, the counting statistics resolution is limited by the RMS width 
in the sample distribution. For our detector array this width is dominated by fluctuations 
in the number of electrons produced at the photo-cathode of the VPD, which is dominated 



by 7-ray counting statistics when t 
noise RMS width is then given by 



re beam is on. During beam on measurements, the shot 



a Ishot = ^2ql v7£, (12) 

where q is the amount of charge created by the photo cathode per detected gamma-ray, / is 
the average photo-current per detector and fs is the sampling bandwidth, set by the 0.4 ms 
time bin width in the time of flight spectrum 19 



IV. ANALYSIS AND RESULTS 
A. Asymmetry Definition 

For a point target and a detector array with perfect spatial resolution, the measured 7-ray 
angular distribution would be proportional to the differential cross section Y = 1 + A 1 cos 8, 
where 9 is the angle between the neutron polarization and the momentum of the emitted 
photon and A lt uD is the parity-odd up-down (UD) asymmetry. A third term is present 
if a parity-conserving (PC) left-right (LR) asymmetry exists |28[. In that case Y = 1 + 
A^jjd cos 9 + A-f^LR sin 9. However, the relationship between the basic expression for the 7- 
ray yield and the measured asymmetry is complicated by a number of small neutron energy 



1(3 



dependent effects. A separate asymmetry is calculated for each detector pair, as defined in 
FigH 




FIG. 4: A ring of detectors and one up-down pair, as seen with beam direction into the page. B 
is the magnetic holding field defining the direction of the neutron polarization. 



The physics asymmetry for a given detector pair (p), spin sequence (j), and neutron time 
of flight (tj) is given by 



+ (Gud(U)) + {A&iU) +PA{% b (t i )) (G LR (U)) 

= {Aj>UU) - Av g A f {U) - Aj mse ) 



(13) 



Here, Ap£ w (ti) is the measured asymmetry. The background asymmetries (A 3 ^ b , A 3 £ P R b ) 
and the relative signal level (/3) must be measured in auxiliary measurements. A p g is the 
gain asymmetry between the detector pair and Af(ti) is the asymmetry from pulse to pulse 
beam fluctuations. The neutron energy and detection efficiency weighted spatial average 
detector cosine (up-down asymmetry) with respect to the (vertical) neutron polarization 
is given by (Gud{U)) — cos(9), while the detector sine (left-right asymmetry) is given by 
{GhRiti)) — sin(8). These detector-target geometry corrections have been modeled for 
each target geometry. Also included are the correction factors due to the neutron beam 
polarization (P n (ij)), the spin flip efficiency (A s ^(tj)) and the neutron depolarization in the 
target (A dep (ti)). 
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The measured (raw) asymmetry (AP r ^ w ) for each detector pair and neutron energy can 
be extracted in the usual way, by forming a ratio of differences between cross-sections to 
their sum. However, to suppress first and second order detector gain drifts j3] the raw 
asymmetries were formed for all valid sequences of 8 macro pulses with the correct neutron 
spin state pattern as shown in eqn. (|14j) . 

A jp _ EU^'fc) ~ D ^ ~ ^s =l (Us(U) ~ Ds(U)) 
raw[ l) Es=,(Us(U) + D.(ti)) + Es =l (Us(U) + D.(U)) ' 1 j 

Here the sum is over all four signals with the corresponding spin state in a spin sequence 
for the up (U) and down (D) detector in a pair. A so-called valid 8 step sequence of spin 
states is defined as (T-lltitt!)- Asymmetries were measured for 55 different neutron energies 
between approximately 2 and 16 meV, with a resolution of ~ 0.2 to 1.0 meV per time bin, 
respectively. 

It is important to realize that signal fluctuations that are not correlated with the switching 
of the neutron polarization direction, such as beam and detector gain fluctuations, will 
average out and don't contribute to the asymmetry. It is, however, essential that these 
signals have an RMS width that is small compared to the RMS width in the asymmetries of 
interest (driven by counting statistics) so that they do not reduce the statistical significance 
of the result and are averaged to zero quickly compared to the time it takes to measure the 
asymmetry to the desired accuracy. Possible false asymmetries due to spin-state correlated 
electronic pickup (additive asymmetry) and possible magnetic field induced gain changes 
(multiplicative asymmetry) in the detector VPDs have previously been measured and are 
consistent with zero to within 5 x 10 -9 jl^ . 

The detector pair physics asymmetries as represented by eqn. ED can then be combined in 
error weighted averages over the neutron time of flight spectrum to form a single asymmetry 
for each detector pair in the array, for a single 8-step sequence of beam pulses. If beam 
intensity levels are sufficiently stable over the measurement time these sequence asymmetries 
can be histogrammed for each pair. Typical run lengths were ~ 8.3 minutes and included 
10000 beam pulses or 1250 8-step sequences and the asymmetry measurements performed 
usually extended over several hundred runs. 
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B. Results 



The known parity-odd gamma asymmetry in CCI4 was used to verify that a nonzero 
asymmetry can be measured with our apparatus. The CCI4 asymmetry was also used to 
verify the geometrical dependence of the pair asymmetries. For this purpose all 24 pair 
asymmetries, extracted from the histogrammed 8-step sequence asymmetries over all data 
obtained with that target, were multiplied by their mean geometry factors and plotted versus 
their corresponding mean angle. The resulting graph is shown in Fig. The fit function 
used to extract the total array asymmetry is Ajj£, cos 9 + Alrsyo. 9. 

■ x10 6 



x1(T 




-0.05b 



2 4 6 8 10 12 14 16 18 20 22 
Pair Angle + ring*2rc (rad) 




-0.15. 



2 4 6 8 10 12 14 16 18 20 22 

Pair Angle + 2ii*ring (radians) 



FIG. 5: Left: CCI4 asymmetries for each pair, plotted versus angle of the first detector in the pair 
w.r.t the vertical. The total array asymmetry is extracted from the fit. Right: Noise asymmetries. 

In general, the up-down and left-right asymmetries must be extracted using the fit de- 
scribed above. Higher order corrections to the fitting function used here (parity violating 
or not) are introduced by higher partial waves in the expansion of the initial and final two 
nucleon states representing more complicated scalar combinations between the neutron spin 
s n and outgoing 7-ray momentum direction k 7 . For the up-down asymmetry the angular 
distribution is obtained from initial and final two-nucleon states with components up to the 
P-waves producing the s n • k 7 correlation. The left-right asymmetry originates from the 
s n • (k 7 x k n ) correlation. Parity violating corrections from higher partial waves are negli- 
gible because they represent a second order perturbation proportional to the weak coupling 
squared. The results of the asymmetry measurements are summarized in table IIVI Note 
that beam asymmetries are only produced if there are pulse to pulse fluctuations in the num- 
ber of neutrons and only in combination with a difference in gain between a given detector 
pair. Neither beam fluctuations nor detector gain differences are correlated with the neutron 



19 



spin and therefore the Beam- Gain asymmetry does not contain any up-down or left-right 
dependence. Due to the sum over the eight step sequence, the Beam- Gain asymmetry is 
zero and its root mean square width is determined by the size of beam fluctuations. The 
additive and multiplicative noise asymmetries in II VI are measured without a light signal 
from the detectors (electronic noise only) and with a light signal from LEDs embedded in 
the detectors respectively. The large RMS width for the multiplicative noise asymmetry is 
a result of larger fluctuations with LEDs [3]. 



C. Errors 

The final statistical errors stated in table IIVI are taken from the distribution of se- 
quence values (J^/N = (E(A^) — E(A 1 ) 2 )/N, with N histogrammed 8-step sequence 
asymmetries. Any non-random effect such as those introduced by the correction factors 
\(G(ti),)\, A dep (ti), P n (ti), Asfl (tj) are treated as systematic errors. These enter as 

— V(f) 2+ te) 2+ (tr+te) 2 



Asymmetries and RMS width 





Up-Down 


Left-Right 


RMS width 
(ty P .) 


Al 


(-0.02 ±3) x 10" 7 


(-2 ±3) x 10" 7 


1.2 x 10" 3 


CC1 4 


(-19 ±2) x 10~ 6 


(-1±2) x 10~ 6 


1.0 x 10~ 3 


B 4 C 


(-1±2) x 10~ 6 


(-5 ±3) x 10" 6 


0.7 x 10~ 3 


Cu 


(-1±3) x 10~ 6 


(0.3 ± 3) x 10~ 6 


1.0 x 10~ 3 


In 


(-3 ±2) x 10~ 6 


(3 ±3) x 10~ 6 


0.4 x 10~ 3 


Noise (add.) 


(2 ± 5) x 10^ 9 


(-7 ±5) x 10~ 9 


2.0 x 10" 6 


Noise (mult.) 


(3 ± 7) x 10~ 9 


(-9 ±7) x 10- 9 


0.2 x 10~ 3 


Beam-Gain 


N/A 


N/A 


1.0 x 10" 5 



TABLE IV: Up-Down and Left-Right asymmetries for the target materials. Stated errors are sta- 
tistical only. The RMS widths are taken from histograms with single 8-step sequence asymmetries 
for a detector pair as individual entries. 
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and are added in quadrature with the statistical error. 

The errors on the beam polarization and spin flip efficiency were calculated to be 4% and 
10% respectively. The error on the geometry factor is estimated to be less than 1% from 
variations observed in the values when varying the step size in the Monte Carlo, simulating 
7-ray interaction in the detectors. The error on the spin flip scattering is estimated to be 
on the order of a few percent. Since the systematic errors are scaled by the asymmetry, their 
contribution to the overall error on the asymmetry is negligible compared to the statistical 
error, except for the case of the CCI4 target, which has a large non-zero asymmetry. For 
CCI4 , the systematic error is ~ 2.3 x 10 -6 . So the CCI4 Up-Down physics asymmetry 
and its total error is (—19 ± 3) x 10~ 6 . A prev ious measurement of this asymmetry by this 
collaboration found (—29.1 ± 6.7) x 10 -6 [3(|. M. Avenier and collaborators found an 
Up-Down asymmetry for 35 C1 of (—21.2 ± 1.7) x 10~ 6 , while V.A. Vesna and collaborators 
found (-27.8 ± 4.9) x 1(T 6 U| (see also 0) 



V. CONCLUSION 



The NPDGamma collaboration has searched for 7-ray asymmetries from polarized slow 
neutron capture on 27 A1, Cu, 115 In and B4C. The asymmetry measurements for these tar- 
gets were consistent with zero at the few 10~ 7 level for 27 A1 and at the few 10 -6 level for Cu 
and 115 In . All asymmetries are consistent with zero within errors. The 35 C1 asymmetries 
obtained from the CCI4 measurements are consistent with results from previous measure- 
ments. A statistical model, in combination with previous measurements of weak matrix 
elements in compound nuclei, was used to estimate the expected RMS size of the parity 
violating 7-ray asymmetries in 27 Al, Cu, and 115 In . Based on this model it is expected that 
non-zero measured asymmetries will be smaller than the estimated width 68.3% of the time. 
The upper bounds on the measured asymmetries are therefore consistent with the estimates 
obtained from these statistical calculations. Based on the inverse relationship between the 
single particle level spacing and the size of the asymmetry, one would expect a large number 
of very small or essentially zero asymmetries when performing measurements for many larger 
nuclei, but one would also expect to find a small number of nuclei with enhanced asymme- 
tries. We plan to continue measurements in other nuclei in the mass range A > 50 to test 
this hypotheses more precisely and to further investigate the predictions of the statistical 
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approach to parity violation in compound nuclei. 
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